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ABSTRACT 

Context. The Lyman-alpha (Lya) recombination line is a fundamental tool for galaxy evolution studies and modern observational 
cosmology. However, subsequent interpretations are still prone to a number of uncertainties. Besides numerical efforts, empirical data 
are urgently needed for a better understanding of the Lya- escape process. 

Aims. We empirically estimate the Lya escape fraction in a statistically significant sample of galaxies in a redshift range z ~ - 0.3. 
This estimate will constrain interpretations of current high-redshift Lya observations. 

Methods. An optical spectroscopic follow-up of a sub-sample of 24 Lya emitters detected by GALEX at z ~ 0.2 - 0.3, combined 
with a UV-optical sample of local starbursts, both with matched apertures, allow us to quantify the dust extinction through Balmer 
lines, and to estimate the Lya escape fraction from the Ha flux corrected for extinction in the framework of the recombination theory. 
Results. The global escape fraction of Lya radiation spans a wide range of values and f esc (Lya) clearly decreases with increasing 
nebular dust extinction E(B-V). Several objects show /„ c (Lya) greater than /^(continuum), which may be taken as observational 
evidence for a clumpy ISM geometry or for an aspherical ISM. Selection biases and aperture size effects may still prevail between 
z ~ 0.2 - 0.3 Lyman-alpha emitters (LAEs) and local starbursts and may explain the difference observed for f esc (Lya). 
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\Q 1 . Introduction line. The determination of the amount of Lya radiation that es- 
capes from the host galaxy is certainly the most important step 

OS , Considerable progress has been made in the last years in the de- toward understanding how various galaxy properties may distort 

O ■ tection and characterization of distant galaxy populations thanks, the interpretation of Lya observations. Although the order of im- 

^ | in particular, to 8- 10m class telescopes with large field of view portance of these parameters has been extensively discussed (e.g. 

. instruments. In this context, the Lya emission line is of par- IVerhamme et alj|2008t fSchaerer & Verhammdl2008T; lAtek et all 

/\ ' ticular interest, since it rem ains the brightest spectral signature 120091) . empirical evidence based upon a large sample of galaxies 

^ I of rem ote young galaxies dPartridge & Peebles] 1 1967t ISchaeretj is still missing. 

- - 1 120031) . As a result, high redshift galaxies are now being rou- ^ T . .. , , T]1 . , . r , » 

. — , , , ■ , ■ T t, i i m ■ From the International Ultraviolet Explorer (IUE) to the 

tinely detected thanks to the Lyman Break selection and/or emis- TT „ „ , , TTC ™, . . , . 

sion line surveys (e g. iGronwall let alJ 2007; Ouchi et al. 2008; H " bble f SpaCe Telesc °P e (Hf) era < spectroscopic and imaging 

rrrr. 71 \ i" h b , - , . ' — — : — . . observations of nearby star forming galaxies have played a key 

Nilsson et al. 2008, and references therein); this situation is , . .- £ ■ . . c T 

\.. , : — J ' ** . , . , — i , ' , role in identifying the main parameters responsible for the Lya 

likely to improve with upcoming Extremely Large Telescopes . , n . , , .. T 

/ctt\ /u i w uu c ti /iwct-\ a escape in a given galaxy. Recent high resolution Lyff imaging re- 

(ELTs) and the James Webb Space Telescope (JWST). A com- r , , , ° . \ . c t c J . ° . 

,' . , , f • , , , r suits clearly demonstrate the importance of resonant scattering, 

parable survey is now available tor the first time at low redshift ., , ; „ , u- u i m j j 

i no a 15 ffTT 1 i MoaaoK ^ i t t i. tatcv evidenced locally by very high Lya/Ha ratios and an outstand- 

^ ° tT °; 3 t 5 '™»IS3 2008 thanks to the GALEX ing large Lya sc ^ tt / ring halo (e. g Hayes et al. 2007; AtekiFd] 

(Galaxy Evolution Explorer) UV capabilities. ^ |g> S tlin et alJl2^ . However, most results so far have no 

The Lyo' line proves an invaluable tool in a cosmological statistical bearing and are still difficult to generalize, because the 

context and is used in a wide variety of applications. Lya has sample is not only small but consists of specific "hand-picked" 

been used in recent studies to probe early stages of galaxy for- objects. We propose here to improve this situation by using a 

mation, estimate the star formation rate, trace large scale struc- larger sample of 24 Lya emitting galaxies at 0.2 SS z ^ 0.35 

tures, identify potential hosts of population III stars and place found by GALEX. We have carried out a spectroscopic follow- 

constraints on cosmic reionization at z i 6. However, in order up of a southern sub-sample with EFOSC2 on the ESO New 

to ensure a proper interpretation of these very promising Lya- Technology Telescope (NTT). This enables us to analyze how 

oriented studies, one first must establish a robust calibration of the Lya emission is related to many physical properties of galax- 

the many parameters that control the complex transport of this ies. We also re-analyzed UV-optical spectra of 1 1 local star- 
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bursts. For the first time, these large aperture observations allow 
us to determine empirically the Lya escape fraction in a large 
sample of galaxies and to examine its dependence on dust ex- 
tinction, if any. 

2. Observations 

2.1. The GALEX sample 

96 Lya emitting galaxies at z ~ 0.2 - 0.35 were found by 
iDeharveng eTaT] (I2008I) in the far ultraviolet (1350A - 1750A) 
from a GALEX slitless spectroscopic survey. Five fields covering 
a total area of 5.65 deg 2 were used to extract all continuum spec- 
tra with a minimum signal-to-noise ratio (S/N) per resolution el- 
ement of 2 in the FUV. Lya emitters are then visually selected 
on the basis of a potential Lya emission feature, which naturally 
leads to a threshold of EWhya ~ 10A. Data re duction and field 
charac teristics are described in more detail in iDeharveng et al.l 
(I2008I) . 

2.2. Spectroscopic follow-up 

Spectroscopic observations of 24 of the 31 galaxies in the 
Chandra Deep Field South (CDFS) and ELAIS-S1 fields 
were performed with EFOSC2 on the NTT at ESO La Silla. 
Observational conditions were very good with photometric sky 
and sub-arcsec seeing (0.5" - 1"). Two instrumental setups were 
used in long slit mode: (1) a spectrophotometric mode with a 5" 
slit, allowing observations to encompass the whole galaxy (20 
out of 24 objects); and (2) a spectroscopic mode with a 1" slit, 
giving a better spectral resolution enabling us to correct Ha data 
for Nil contamination. Both settings were used in combination 
with Grism #13, covering a large wavelength range in the op- 
tical domain (3690 - 9320A). A binning of 2 x 2 is used and 
corresponds to a plate-scale of 0.24" px _1 and a spectral reso- 
lution of FWHM ~ 12 A (for 1" slit spectra). To avoid second 
order contamination that affects the longer wavelength range, an 
order sorting filter has been mounted to cut off light blue-ward 
of 4200 A. Observational settings and the mean exposure time 
per object are summarized in Table Q] 



Table 1. NTT EFOSC2 observations. 



Mode 


Slit 


Exptime (sec) 


Nb Obj 


Grism 


Spectroscopy 


1" 


3000 


24 


13 


Spectrophotometry 


5" 


1800 


20 


13 



The EFOSC2 spectra were reduced and calibrated using 
standard IRAF routines. The aperture extraction of ID spectra 
was performed through the DOSLIT task. Spectra were flux cal- 
ibrated using a mean sensitivity function determined by obser- 
vations of stan dard s tars (Feigel 10, HILT600, LTT1020, EG21) 
from the lOkd d 19901) catalog. 

2.3. IUE starburst sample 

We have re-analyzed UV-optic al spectra of 1 1 local star- 
burst galaxies, presented in iMcOuade et alJ JT995) and 



IStorchi-Bergmann et ail d 19951) (see also lGiavalisco et alJ[T996l) 
that are distant enough to separate the Lya feature of the galaxy 
from geocoronal Lya emission. In this way, we have com- 
plete control on the Lya emission measurement method. Indeed, 
the definition of a Lya emitter/absorber could be ambiguous 
for P Cygni profiles or an emission blended with absorption. 
Therefore, we consider here only the net Lya flux. Data consist 
of IUEUV (1200-3300 A) spectra combined with ground-based 
optical spectra with a matched aperture (20" x 10"). The line 
measurements have been performed following the same proce- 
dure used for the NTT spectra. We included in our analysis only 
net Lya emitters, i.e. with E Wi^ a > 0. Furtherm ore, for common 
objects between IUE and the I Atek et al.l d2008l) samples, we re- 
tained the latter to minimize aperture size effects (see below). 
We then obtained five IUE objects. 

2.4. Emission line measurements 

All spectra were analyzed using the SPLOT package in IRAF. The 
redshift was measured using the wavelength position of several 
lines and line measurements were performed interactively on 
rest-frame spectra. We confirm with a better accuracy the red- 
shift determinations based on the blind search for Lya features 
in the GALEX survey. 

Fluxes and equivalent widths (EWs) were measured for Ha, 
H/3 and [Nn] 6548, 6584 A. For most spectra, the Ha line (6563 
A) is blended with [Nn] lines, even for the 1 " slit observations. In 
this case, a deblending routine is used within SPLOT to measure 
individual fluxes in each line. Then, the Nn/Ha line ratio is used 
to correct the spectrophotometric observations for [Nn] contam- 
ination. It appears that the dust extinction (cf. Sect. 12.51 ) is sensi- 
tive to the aperture size, since our 1" slit observations, targeting 
the center of the galaxies, led in general to a higher extinction. 
But this does not affect the Nn correction, which is relatively in- 
sensitive to dust extinction. To correct for underlying stellar ab- 
sorption in the Balmer lines, we assumed a constant equivalent 
width of 2 A, typical for star burst galaxies dTresse et al.l l 1996; 
iGonzalez Delgado et alii 19991) . 

To determine uncertainties in the line fluxes, we ran 1000 
Monte Carlo simulations in which random Gaussian noise, based 
on the data noise, is added to a noise-free spectrum. Then, emis- 
sion lines were fitted. The computed MC errors depend essen- 
tially on the S/N quality of spectra. Error propagation is applied 
through the calculation of all the quantities described above and 
the line ratios, extinction etc, computed hereafter. 

Using BPT diagr ams dBaldwin et al l 119811; 
IVeilleux & Osterbrockl 1 19871) . and Chandra X-ray observa- 
tions, we have identified at most three galaxies possibly excited 
by an active galactic nucleus (AGN), which represents up 
to 12.5 % of our sample. Studying similar z ~ 0.3 samples, 
Scarla ta et al. (2009, submitted) find a comparable value (17 %), 
while Finkelst ein et alJ d2009al) claim a much higher fraction 
(around 43 %). We therefore removed these objects from our 
analysis. This diagnostic will be thoroughly addressed in Atek 
et al. (in prep). 

2.5. Extinction 

Reddening along the galaxy line of sight is caused by interstel- 
lar dust extinction. The reddening contribution of our Galaxy 
is negligible for our objects. Then, the extinction coefficient, 
C(H/3), intrinsic to the observed object can be calculated using 
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the Balmer ratio between Ha and H/J: 

f (Ha "> = Rx \Q-C[S(Ha)-S(HB)] qn 

where f(Ha) and f(H/3) are the measured integrated fluxes and R 
is the intrinsic Balmer ratio. We use here a value of R = 2.86, 
assu ming case B rec ombination theory and a temperature of 10 4 
K (Osterbrock 1989). S (Ha) and S (Hp) are determined from the 
Cardelli et al. (1989) extinction law. The colour excess Eb-v is 
then simply computed using Eq. [T]and the relation E(B - V) = 
C/1.47. Due to photometric calibration errors and probably a 
stronger stellar absorption than assumed here, few objects have 
a negative extinction with large uncertainties. This could also be 
the result of enhanced H/3 emission from a reflection nebula. 

3. Lya escape fraction 

In Fig. [1] the Lya flux is plotted against the Ho- one for the 
GALEX objects . While some objects show a Lya/Ha ratio con- 
sistent with the case B value, most of them lie well below this 
line. On the other hand, few objects show a Lya/Ha > 8.7 at a 
significant level. The heavily attenuated values can be explained 
by the high extinction in the UV compared to the optical and to 
resonant scattering of Lya that increases the absorption optical 
depth. However, the observed dispersion and the points above 
the case B line could be the result of other parameters that will 
be discussed below. A convenient way to constrain the dust ex- 
tinction parameter is to determine the escape fraction of Lya as 
a function of E(B-V). 




f(Ha) [erg cm -2 s"'] 

Fig. 1. Lya versus Ha fluxes for the GALEX sample. The dashed 
line represent Lya/Ha=8.7 (case B) and Lya/Ha=l. The local 
objects (not represented here) show higher fluxes and lie in gen- 
eral close to the line 1:1. 

To determine the Lya escape fraction we follow lAtek et alj 
(120081) . The method relies on the fact that Ha emission is not 
prone to complex radiation transport effects but is only affected 
by dust attenuation. Therefore, correcting the observed Ha flux 
for extinction wh ile assuming a case B recombination theory 
dOster brock 1989), one can estimate the intrinsic Lya flux. The 
Lya escape fraction is then given by: 

f esc (Lya) = /(Lya)/(8 .7 x /(Ha) c ), (2) 




-0.2 0.0 0.2 0.4 0.6 0.8 

E(B-V) 



Fig. 2. Lya escape fraction as a function of dust extinction, 
observed in z ~ 0.3 Lya galaxies. The red dashed line rep- 
resents the best fit to our entire sample (GALEX, IUE, and 
Atek et al. samples). The dark dot-dashed line represent the 
best fit determined by Verha mme et alj d2008l) from spectral fit- 
ting of z ~ 3 LBGs. The blue dotted line corresponds to the 
escape fractio n of the cont i nuum attenuated by dust extinc- 
tion using the ICardelli et alJ <t!989h law. Green points are pre- 
dictions for f esc (Lya) using a 3D Lya radiation transfer code 
d Verhamme et al.ll2008l) . See text for details. 

where /(Lya) is the observed flux and /(Ha)c is the extinction- 
corrected Ha flux. 

Figure [2] shows our empirical f esc values as a function of the 
nebular extinction for the GALEX and IUE samples describe d 
above. f esc values for 4 nearby galaxies from I Atek et al.l (|2008), 
with EWi,y„ > A, are also shown. This figure summarizes much 
information with several implications for Lya physics. 

First, we find that f esc is not constant. It spans a wide range of 
values, typically from / esc ~ 0.5 to 100 % in the GALEX sample. 
Second, the Lya escape fraction is clearly sensitive to the dust 
extinction and an anti-correlation is observed. We performed a 
linear least-squares fit to this anti-correlation. The red dashed 
curve is the best solution found corresponding to 

/ MC (Lya) = MLya) E(B-V) . k(Lya) ^ i 2 J + OA. (3) 

The extinction coeffici ent at Lya wavele ngt h goes from 
£(121 6) ~ 9.9 to 12.8 for lCardelli et alj (fT989h and lCalzetti et all 
(2000) laws, respectively. Here, k(Lya) derived from the fit takes 
into account the averaged effects of all processes affecting the 
Lya escape, such as the resonant scattering experienced by Lya 
photons, which increases their mean path and therefore the effec- 
tive dust optical depth, or velocity fields in the gas or the ISM ge- 
ometry, which may ease the escape of Lya. The observed scatter 
around our mean relation is most likely indicative of this multi- 
parameter process. Us ing hydrodynamica l cosmological simula- 
tions of Lya emitters, iDaval et all d2009l) find similar trend be- 
tween f esc (Lya) and extinction. For similar extinctions, the Lya 
escape fraction tends to be higher on average in GALEX Lya 
galaxies than in local objects. In the z ~ objects, f esc never ex- 
ceeds ~ 10%, whereas it covers a wide range of values reaching 
100% at z ~ 0.3. This difference may be the result of different se- 
lection effects. The GALEX galaxies are selected from their Lya 
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emission amongst spectra taken from a blind search, whereas 
the local objects are from a sample of specific, optically-selected 
starburst galaxies. 

Third, several objects show / eje (Lya) greater than the es- 
cape fraction expected for the continuum near Lya (blue dotted 
line) as given by / esc (cont) = 10~ 04 k(12l6> E{B - V \ where the most 
favourable value £(1216) ~ 9.9 is adopted from the Cardelli et 
al. law. Although Lya/Ha ratios exceeding th e theoretical valu e 
have already been found in local starbursts (lAtek et al.l 12008), 
they only occur locally, in spatially resolved objects, where this 
easily can be explained by a local Lya "excess" due to scat- 
tering. Here, in some objects, it is the "global" Lya escape 
fraction determined from the integrated spectra that is found 
to be higher than expected from the most favourable (i.e. flat- 
test) attenuation law. These objects may be obs ervational evi- 
dence for a multi-phase configuration of the ISM (Neufeld 199U 
lHansen & Peng Ohll2006l iFinkelstein et alj|2009bl) . where dust 
is primarily distributed in cold neutral clouds with an ionized 
inter-cloud medium. By reflecting on the cloud surface, Lya 
photons will be easily transmitted through the ionized medium. 
Alternatively, Scarlata et al. (2009, submitted) have advocated 
for a clumpy dust distribution scenario able to reproduce the ob- 
served Lya, Ha and H/3 intensities without the need for a two- 
phase model responsible for different paths for Lya and Ha. 

Approximately 2/6 of these objects also show a relatively 
high Lya equivalent width (EWi_, ya ~ 100-150 A), as may be ex- 
pected for a clumpy ISM. Furthermore, / esc (Lya) higher than 
the UV continuum could also be due to orientation effects in ob- 
jects with an aspherical ISM, e.g. in conical outflows, into which 
Lya would be "channeled" more effectively than continuum ra- 
diation. The Lya escape fractions of objects below the attenu- 
ation curve of the continuum can be quantitatively reproduced 
by models using a homogeneous ISM. However, this does not 
exclude that clumping may also play a role in these objects. 

We performed an extensive grid of 3D Lya radiation trans- 
fer simulations in homogeneous, spherically expanding shells 
(Hayes et al. in prep), around a central emitting source of UV 
continuum plu s the Lya line, using an updated version of the 
MCLya code dVerhamme et al.ll2006l) . We overplot in Fig.|2]the 
predictions of f esc as a function of E(B-V) by using all possible 
values of the remaining parameters affecting f esc : the expansion 
velocity of the neutral gas, Hi column density, Doppler param- 
ete r b and FW HM(Lya) of the input emission line (for details, 
see Atek et al.ll2009h . The model E(B-V) is computed from the 
predicted UV continuum attenuation, assuming the same extinc- 
tion law as above. The model grid covers the bulk of the ob- 
served variations in f esc and E(B-V). The homogeneous, spheri- 
cal shell models may in principle be able to explain the majority 
of the objects, although tailored models including all observa- 
tional constraints are needed to confirm this. However, the f esc 
values above the continuum attenuation curve (blue dotted line 
in Fig. |2]i observed in 6 objects cannot be explained with these 
models, since Lya photons cannot be less attenuated than the 
continuum in a homogeneous ISM. New radiation transfer com- 
putations in clumpy media are underway to examine this inter- 
esting behavior. 

4. Discussion 

We have presented here an estimation of the mean Lya escape 
fraction as a function of the extinction and how different param- 
eters can alter this simple relation. We have carefully chosen 
our aperture (5" slit) in optical spectroscopy in order to obtain 
Balmer fluxes comparable to Lya ones obtained by the GALEX 



grism. This allows us to estimate the extinction and the escape 
fraction for the whole galaxy. 

To keep the same consistency between the different samples, 
we decided to minimize selection effects by retaining only the 
net Lya emitters in both the IUE and the imaging samples. We 
recall that the retained IUE large aperture (20" x 10") UV ob- 
servations and ground-based optical spectroscopy of local star- 
bursts are all aperture matched. Furthermore, for Lya imaging 
objects, the HST aperture is large enough to encompass the en- 
tire Lya emission region. However, while the slitless mode of 
GALEX enables us to recover the diffuse Lya emission, this is 
not necessarily the case for IUE observations. The size of the 
large aperture may remain insufficient for some nearby objects 
in order to encompass the scattered photons across a large area 
of the galaxy, as usually indicated by the large extent of the Hi 
gas. Therefore, / csc (Lya) could be slightly underestimated. On 
the contrary, as mentioned earlier, the z ~ 0.2 - 0.3 objects are 
selected on the grounds of their Lya strength. This will likely fa- 
vor the high Lya escape fractions. One should then keep in mind 
that these opposite effects contribute to stretch the deviation 
from our best fit of this compilation of data points. Similar Lya 
equivalent width criteria are commonly applied to select high-z 
LAEs, implying most likely relatively high escape fractions and 
a less severe discrepancy between Lya and non-resonant radi- 
ation for this class of objects. A blind search of Lya emitters 
would then find lower f esc (Lya) than the GALEX objects. This 
is already found in our current double blind survey in Lya and 
Ha at z ~ 2.2, where an average escape fraction of ~ 4.5 % is 
measured (Hayes et al. in prep). 

In contrast with other emission lines, the dust extinction is 
only one of several parameters governing / fJC (Lya), and the ex- 
tent of the dispersion around the fit is a good illustration. Given 
the importance of an accurate estimation of / fJC (Lya), one needs 
to quantify the kinematical effects by measuring the cold ISM 
velocity in these objects; an additional step toward a precise cal- 
ibration of / esc (Lya). A more detailed investigation of the phys- 
ical properties and SED modeling of our sample will be carried 
out in subsequent publications (Atek et al. in prep). 
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